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ABSTRACT 

Areas  of  convective  precipitation  were  delineated 
and  averaged  convective  precipitation  amounts  determined 
within  these  areas.   A  statistical,  step-wise,  screening, 
linear  regression  procedure  was  used  to  correlate  certain 
large  scale  parameters  with  the  convective  precipitation 
amounts. 

Of  the  parameters  considered  in  the  two  cases 
investigated,  the  best  predictor  of  convective  precipita- 
tion was  either  the  terrain-induced  vertical  motion  at 
the  top  of  the  boundary  layer,  or  the  product  of  the 
latter  and  the  specific  humidity  at  the  same  level. 
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1.   INTRODUCTION 

The  failure  of  numerical  models  to  account  for 
convective  precipitation  is  well  known  (see,  e.g., 
Danard  (1966  a,  1966  b)).   Such  a  shortcoming  is  expected 
since  the  equations  used  for  numerical  prediction  are  not 
applicable  to  sub-synoptic  scale  motion.   Convective 
precipitation  is  highly  variable  in  both  time  and  space. 
This  would  prevent  such  precipitation  from  being  accu- 
rately accounted  for  even  with  a  much  denser  observational 
network  than  the  one  available  in  the  United  States  today. 

This  investigation  is  an  attempt  to  delineate  the 
areas  of  convective  precipitation  by  the  use  of  semi- 
empirical  methods.   Convective  precipitation  will  be 
interpreted  here  as  the  average  amount  associated  with 
sub-synoptic  scale  systems  over  a  square  area  whose  side 
is  equal  to  the  distance  between  adjacent  grid  points 
used  in  numerical  prediction  (381  kilometers  at  60  degrees 
north).   The  amount  defined  is  assumed  to  be  determined 
by  the  large-scale  properties  of  the  atmosphere. 

Convective  activity  is  confined  to  areas  of  unstable, 
or  conditionally  unstable,  lapse  rates.   That  is,  the  lapse 
rate  should  be  steeper  than  the  moist  adiabatic  lapse 
rate.   Kuo  (1965)  describes  a  convective  cloud  model  which 
will  "bring  surface  air  to  all  levels  up  to  a  great 
height  so  that  inside  the  cloud  the  vertical  distribution 
of  temperature  and  mixing  ratio  are  those  of  the  moist 


adiabat  through  the  appropriate  condensation  level." 
This  suggests  that  the  rate  of  production  of  cloud  and 
total  water  content  would  depend  on  the  specific  humidity, 
or  mixing  ratio,  and  the  vertical  velocity  at  some  lower 
boundary.   Ooyama  (1963)  assumed  that  the  rate  of 
precipitation  was  proportional  to  the  flux  of  water 
vapor  through  the  top  of  the  frictional  boundary  layer. 
The  total  gain  of  water  vapor  in  an  air  column,  assuming 
no  flux  through  the  top  of  the  column,  may  be  represented 
by 

(ecg.,  Krishnamurti  (1968)),  where  I  is  the  water  gained, 
q  is  the  specific  humidity,  and  p  is  pressure.   The  sub- 
scripts b  and  t  indicate  the  bottom  and  top  of  the  column, 
respectively.   Here  the  first  term  represents  the  flux  of 
water  vapor  through  the  sides  of  a  column  of  unit  area, 
and  the  second  term  represents  the  flux  through  its  base, 
which  is  taken  here  as  the  top  of  the  planetary  boundary 
layer.   Since  specific  humidity  decreases  rapidly  with 
height,  the  second  term  is  frequently  larger  than  the 
first  if  the  level  p,  is  sufficiently  high. 
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2 .   PROCEDURE 

For  the  two  case  studies  included  here,  Danard 
(1966b)  has  made  numerical  predictions  of  gridpoint 
precipitation  within  an  area  shown  in  Figures  5-8* 
Predicted  amounts  obtained  from  the  first  12  hours  of 
these  integrations  were  assumed  to  represent  precipitation 
effects  of  large-scale  motions.   Observed  precipitation 
amounts  were  obtained  by  averaging  all  individual  values 
reported  by  synoptic  stations  over  a  square  centered 
on  the  grid  point  and  having  a  side  equal  to  the  distance 
between  adjacent  grid  points  (approximately  381  kilometers) 
Convective  precipitation  amounts  used  are  the  difference 
between  the  averaged  observed  amounts  and  the  amounts 
predicted  by  Danard,   Only  those  points  where  the  lapse 
rate  was  close  to,  or  greater  than,  the  pseudo-adiabatic 
value  (see  below)  were  used.   The  maximum  12 -hour 
precipitation  amount  predicted  by  the  numerical  scheme 
in  these  areas  was  only  0.4  millimeters  in  either  case 
study.   Thus,  the  areas  of  convective  precipitation  were 
clearly  delineated  from  those  of  large-scale  precipitation. 

A  statistical,  step-wise,  linear  regression  procedure 
was  used  to  correlate  convective  precipitation  with  the 
following  parameters:  (l)  the  difference  between  the  pseudo- 
adiabatic  and  the  computed  lapse  rates  (Q~0,    (2)  the 
vertical  velocity  (oj)  at  the  top  of  the  boundary  layer, 
(3)  specific  humidity  (q)  at  the  top  of  the  boundary 
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layer,  (4)  the  product  of  (2)  and  (3),  and  (5)  dew  point 
depression  (T-T, )  in  the  lower  troposphere.   The  program 
used  is  available  in  the  Naval  Postgraduate  School 
computer  library,  where  it  is  listed  as  BMD02R.   A 
sequence  of  multiple  linear  regression  equations  is 
computed  in  a  step-wise  manner.   At  each  step  the 
variable  is  added  which  makes  the  greatest  reduction 
in  the  error  sum  of  squares. 

For  the  two  cases  considered,  average  temperatures 
and  dew-point  temperatures  were  computed  for  the  85O, 
700,  and  500  mb  levels.   The  reported  values  at  the 
beginning  and  end  of  the  period  were  used  to  compute 
time  averages.   These  averaged  temperatures  were  used 
in  the  computations  . 

All  averaged  parameters  were  computed  from  values 
reported  by  the  stations.   The  resulting  fields  were  then 
analyzed,  and  the  values  for  the  grid  points  used  were 
taken  from  the  analysis. 

As  used  here,  the  lapse  rate  is  defined  to  be  the 
difference  between  the  temperatures  at  85O  and  500  mb . 
The  pseudo-adiabatic  lapse  rate  between  85O  and  500  mb 
was  taken  from  an  Arowagram  as  a  function  of  the  tempera- 
ture at  700  mb „   Convective  precipitation  was  considered 
possible  where  the  difference  between  the  lapse  rates 
(rw~r)  was  less  than  2.0  C/(350  mb ) .   The  areas  delineated 
by  this  criterion  included  all  regions  where  large  amounts 
of  convective  precipitation  was  observed  in  both  cases. 
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The  dew  point  depression  used  was  averaged  between 
85O  and  700  mb.   Here  the  500  mb  level  was  omitted 
because,  in  both  cases  many  stations  reported  motor- 
boating  at  that  level.   In  the  cases  where  motor- 
boating  was  reported  at  the  lower  levels,,  the  minimum 
dew  point  depression  for  which  a  radiosonde  normally 
reports  motor-boating  was  used. 

Danard  (1966a)  expressed  the  vertical  velocity  at 
the  top  of  the  boundary  layer  as 

O)   z:  oOr    +    U)f  D  ( 2  ) 

where  00r  and  6*J»f  are  the  vertical  velocity  induced  by 
orographic  and  eddy-stress  influences,  respectively. 
In  the  area  considered  in  this  study  (see  Figure  l) 
most  of  the  terrain  heights  are  less  than  400  meters 
above  sea  level.  Therefore,  the  surface  geostrophic 
wind  was  taken  to  be  that  at  1000  mb .  Thus,  u>t  is  well 
approximated  by 

where  p,  is  the  standard  atmosphere  pressure  at  the 
elevation  of  the  smoothed  terrain  provided  by  Berkofsky 

and  Bertoni  (1955)*  and  (  V"  )1000  is  the  Seos1:roPhic 
wind  at  1000  mb . 
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The  frictionally  induced  vertical  velocity  was 
computed  by 

uif  a  %*?  CK'Vx  cd  Vj  v9 )      •       (4) 

Here  V  is  a  representative  geostrophic  wind  in  the 

o 

planetary  boundary  layer  and  C,  is  the  drag  coefficient 
(Cressman  (i960)).   In  the  cases  considered  here  (  V  )iqqq 
was  used.   However,  over  higher  terrain,  or  in  cases 
where  (  V  )i0qq  cannot  be  easily  computed,  (  V  )oc-0 
might  be  more  useful. 

Specific  humidity  at  the  top  of  the  boundary  layer 
may  be  expressed  as 

7  -  p-  0.3735         5  (5) 

where  e  is  the  vapor  pressure  and  p  is  the  pressure  in 
mib  at  the  top  of  the  boundary  layer,  assumed  to  be  85O 
mb  here.   Vapor  pressures  used  were  the  saturation 
values  (Smithsonian  Institution  (1951))  corresponding 
to  the  dew  point  temperatures  at  850  mb . 
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3.   RESULTS  OF  CASE  STUDIES 

Two  synoptic  situations  were  investigated:  1200  GMT 
11  February  to  0000  GMT  12  February,  1965*  and  1200  GMT 
16  March  to  0000  GMT  17  March,  1965.   These  intervals  will 
be  referred  to  as  Cases  A  and  B,  respectively.   The 
surface  synoptic  situation  for  the  beginning  and  end  of 
the  two  12  hour  periods  may  be  seen  in  Figures  1-4. 
Both  cases  were  characterized  by  deepening  cyclones 
moving  from  the  southwest.   These  storms  produced  some 
of  the  worst  weather  of  the  1965  winter  in  the  central 
United  States. 

The  number  of  grid  points  used  in  Cases  A  and  B 
(i.e.,  the  sample  sizes)  was  16  and  17,  respectively. 
Tables  1  and  2  give  the  simple  correlation  coefficients 
between  all  the  variables  used  for  Cases  A  and  B, 
respectively.   In  both  cases  the  correlation  between  CO 
and  o)q  was  very  high.   In  Case  A,  o)q  had  the  highest 
simple  correlation  with  convective  precipitation  (P). 
For  Case  B,  the  correlation  between  cO  and  P  was  slightly 
higher  than  that  between  u)q  and  P.   The  second  highest 
simple  correlation  for  Case  A  was  between  (T-T,)  and  P. 
However,  there  was  little  or  no  correlation  between  these 
variables  in  Case  B. 

The  results  of  the  step -wise  regression  analysis  may 
be  seen  in  Tables  3  and  4  for  Cases  A  and  B,  respectively. 
In  the  step-wise  regression  procedure,  the  first  variable 
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entered  is  always  the  one  with  the  highest  simple 
correlation  with  P,  unless  the  "deletion  option"  of 
the  step-wise  regression  program  is  invoked.   In  Case 
A,  after  o)q  entered  the  regression  equation,  u)  was 
essentially  redundant  and  did  not  enter  the  equation 
at  all.   For  the  same  reason  a)q   did  not  enter  the 
equation  in  Case  B,  since  a)  entered  first.   In  both 
cases  there  was  little  improvement  in  the  multiple 
correlation  between  the  predicted  value  and  P  after 
the  first  variable  entered  the  regression  equation. 

An  attempt  was  made  to  determine  if  a  time  delay 
between  the  flux  of  moisture  through  the  boundary 
layer  and  resulting  precipitation  could  be  detected. 
To  accomplish  this,  the  regression  procedure  was 
repeated  using  variables  computed  for  the  beginning, 
rather  than  for  the  average,  of  the  time  period.   The 
resulting  simple  correlation  coefficients  are  given  in 
Table  6  and  7.   In  Case  A,  there  was  noticeable  improve- 
ment in  the  correlation  between  both  o)q  and  P,  and  lP 
and  P.   However  for  Case  B  there  was  little  change  in 
these  correlations. 

The  step-wise  regression  program  used  permits  the 
variable  oJq   to  enter  the  regression  equation  first  for 
Case  B  if  the  variable  lO  is  deleted.   In  this  way  it  is 
possible  to  obtain  the  variable  o>q  as  the  primary 
predictor  for  both  Cases  A  and  B.   In  a  similar  manner,  a) 
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has  been  selected  as  the  primary  predictor  for  Case  A. 
The  resulting  equations,  with  only  one  variable  entered, 
may  be  seen  in  Table  5.   Due  to  the  high  simple  corre- 
lation between  u)  and  u)q  there  is  little  loss  in  explained 
variance  for  either  case,  when  ca)  or  a)q  is  employed  as 
the  primary  predictor. 

In  order  to  improve  the  prediction  equations  and 
to  ensure  that  P2,0,  predictions  were  made  for  both 
cases  using  the  equation 

P  =  A0  +  Ai(^)  ,  for    a>q<-A0/A1  ,^  * 

=   O  3  for    ooq^-Ao/^t 

The  coefficients  A  and  A.,  were  determined  by  computing 


the  quanity 

S2-  I(Pob-P)2        ,  (7) 

where  the  summation  extends  over  all  the  points  in  the 
sample.   P  ,  is  the  observed  convective  precipitation  at 

a  gridpoint  and  P  is  the  amount  given  by  Equation  (6). 

2 
S  is  simply  the  error  sum  of  squares.  Various  values 

p 

of  A  and  A.,  were  tried  until  the  minimum  in  S  was 

found.   This  procedure  was  repeated  using  a)  as  the 
predictor.   The  improved  prediction  equations  are  shown 
in  Table  9.   Comparing  the  last  columns  of  Tables  5  and 
9,    it  is  seen  that  the  above  procedure  results  in  a  slight 
reduction  in  error  for  Case  B.   However,  for  Case  A  the 
improvement  is  negligible. 
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Predictions  were  made  for  both  cases  using  the 
equation 

p  =  -0.45  -  i.ZO(u)q)   3  for  a)Q<-  0.375      ,    x 

(8) 
=  O  jfor   o)^2-  0.37 5  . 

Equation  (8)  was  obtained  by  averaging  the  coefficients 
given  in  the  first  two  lines  of  Table  9.   The  results 
may  be  compared  with  the  observed  convective  precipitation 
in  Figures  5-8. 


18 


Figure   1 


Figure   2 
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TABLE  I 


Matrix  showing  the  simple  correlation 
coefficients  between  the  variables  for 
Case  A,  P  is  convective  precipitation. 


Tw-r     fld        q         u)q     T-T 


d 


P      1.00 

r„-r 

(A) 

q 

0)q 


0.28 

-0.43 

0.12 

-0.51 

-0.47 

1.00 

-0.34 

-0.13 

-0.35 

-0.44 

1.00 

-0.10 

0.95 

0.55 

1.00 

-0.28 
1.00 

0.01 
0.53 

TABLE  II 
Same  as  Table  I   for  Case  B 


r„-  r  u)  q  cjq  T-T 


d 


P  1.00 

rw-r 

q 

uOq 


-0.25 

-0.67 

-0.24 

-0.65 

-0.08 

1.00 

0.01 

-0.65 

-0.01 

-0.34 

1.00 

0.45 

0.98 

0.31 

1.00 

0.48 
1.00 

0.36 
0.28 
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TABLE  V 


Regression  equations  using  only 
the  first  variable  entered. 


Case 

Variable 
entered 

R 

Regres 
coeffi 

Ao 

sion 
cients 

Al 

Root  mean 
square  error 
of  estimate 

A 

a>q 

0.51 

0.57 

-1.66 

4.65 

B 

uJq 

0.65 

-0.01 

-0.32 

1.32 

A 

o) 

0.43 

1.07 

-9.20 

4.89 

B 

u) 

0.67 

-0.08 

-1.37 

1.28 

TABLE  VI 


Matrix  showing  the  simple  correlation 
coefficients  between  the  variables  for 
Case  A.   The  primed  variables  were  computed 
for  the  beginning  of  the  time  period. 


c-r 


a) 


a>V 


t-t' 

d 


r„-r 


toV 


1.00 


0.28 

-0.65 

0.04 

-0.64 

-0.10 

1.00 

-0.24 

-0.39 

-0.14 

-0.30 

1.00 

-0.40 

0.99 

0.36 

1.00 

-0.42 
1.00 

-0.62 
0.33 
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TABLE  VII 


Same  as  Table  VI  for  Case  B 


r„-r 


u>' 


a>'q' 


T'< 


1.00 


rw-r 


u>'q' 


0.25 

-0.62 

-0.11 

-0.66 

0.37 

1.00 

0.06 

-0.64 

0.04 

-0.04 

1.00 

0.44 

0.93 

-0.32 

1.00 

0.41 
1.00 

0.01 
-0.26 

TABLE  VIII 


Same  as  Table  V  for  variables 
computed  for  the  beginning  of 
the  time  period. 


Case 

Variable 
entered 

R 

Regression 
coefficients 

Ao      Al 

Root  mean 
square  error 
of  estimate 

A 

cd'q' 

0.64 

1.71 

-2.08 

4.16 

B 

a)'q' 

0.66 

0.21 

-0.72 

1.31 

A 

CO' 

O.65 

2.22 

-13.77 

4.14 

B 

oO' 

0.62 

0.07 

-2.59 

1.36 
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TABLE  IX 


Improved  regression  equations  usin^ 
only  the  first  variable  entered. 


Regressi 

on 

Root  mean 

Variable 

coeffici 

ents 

square  error 

Case 

entered 

Ao 

Al 

of  estimate 

A 

uJq 

0.20 

-1.80 

4.63 

B 

u)  q 

-1.10 

-0,50 

1.19 

A 

u) 

1.00 

-0.35 

4.88 

B 

u) 

-0.80 

-1.85 

1.20 
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4.   CONCLUSIONS  AND  REMARKS 

Probably  the  most  that  could  be  claimed  for  the 
results  of  this  study  is  that  it  may  have  some 
applicability  over  the  central  and  eastern  United  States 
for  the  winter  months  during  the  period  1200  to  0000  GMT 
However  there  is  indication  that  the  product  of  vertical 
velocity  and  specific  humidity  at  the  top  of  the 
boundary  layer  can  be  used  to  predict  convective 
precipitation . 

While  this  study  is  diagnostic  in  nature,  it  uses 
data  that  can  be  predicted  numerically.   Consequently, 
the  technique  could  be  included  in  a  prognostic  scheme 
with  some  expectation  of  success.   The  technique  of 
using  lapse  rates,  as  described  here,  to  determine  areas 
of  possible  convective  activity  could  be  applied  to 
local  forecast  areas  using  forecast  temperature  fields. 

For  a  synoptic  situation  to  be  used  in  this  type 
of  a  study,  there  must  be  a  clear  delineation  between 
large-scale  and  convective  precipitation.   To  find  out 
if  a  situation  meets  this  requirement,  a  prediction  of 
precipitation  amounts  due  to  large-scale  processes  must 
be  made  and  compared  to  the  reported  precipitation.   The 
time  required  for  this  procedure  limited  this  study  to 
the  two  cases  considered  here.   Furthermore,  the  small 
size  of  the  samples  virtually  precludes  a  statistical 
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proof  of  the  significance  of  the  results.   However, 
the  findings  of  this  study  are  sufficiently  encouraging 
to  justify  further  work. 
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